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We address the study of controllability of a closed quantum system whose dynamical Lie algebra is 
generated by adjacency matrices of graphs. We characterize a large family of graphs that renders a 
system controllable. The key property is a novel graph-theoretic feature consisting of a particularly 
disordered cycle structure. Disregarding efficiency of control functions, but choosing subfamilies of 
sparse graphs, the results translate into continuous-time quantum walks for universal computation. 
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I. INTRODUCTION 

The study of classical control theory is ubiquitous 
across engineering disciplines. The concept of controlla- 
bility at the quantum level is a fundamental notion which 
expresses the ability of implementing any dynamics in a 
given quantum mechanical set-up (see |2| for an intro- 
ductory monograph). From the practical point of view, 
the successful realization of quantum devices for a vari- 
ety of information processing tasks strongly depends on 
the ability of manipulating systems with sufficient free- 
dom. The design of protocols to control closed quantum 
systems mainly deals with schemes for efficient control- 
lability by acting on subspaces ||. 

A large number of systems have been shown to exhibit 
characteristics that allow controllability. For example, 
almost any pair of Hamiltonians that can be coherently 
applied to a finite-dimensional quantum system renders 
it controllable, and almost any quantum logic gate is uni- 
versal [|J. However, controllability of large systems does 
not directly give an efficient implementation of control 
functions, and therefore processes like universal quantum 
computation. Moreover, control criteria are generally not 
computable for large systems and are not immediately 
scalable with the system size @. 

Here we consider controllability in relation to a class 
of dynamics that can be interpreted as Schrddinger evo- 
lutions of a particle hopping between the vertices of a 
graph. The corresponding Hamiltonians are matrices 
with nonzero entries only where transitions are permit- 
ted by the graph structure. This class naturally em- 
braces continuous-time quantum walks (for short, QWs) 
Q, and the evolution of a single (or multiple Q) exci- 
tation subspace for systems of spin-half particles with 
various kind of interaction (e.g., XY, XYZ, etc.) [l7j . 
QWs have been shown to build gates by scattering off a 
set of small graphs attached to wires representing basis 
states @. QWs give examples of matrices for sufficient 
control, when we can arbitrarily modify edge- weights (or, 
equivalently, strength of couplings) jlO ]. 

During our discussion, an n-level system is controllable 
by a given set of Hamiltonians (possibly acting on specific 
subsystems only) if every element of the unitary group 
U(n) can be approximated by the matrices of the sub- 



group obtained via the dynamical Lie algebra of the set. 
This general definition is useful to isolate the main differ- 
ence between the notions of controllability and universal- 
ity. For instance, global evolution of a spin system may 
require complex protocols to implement 2-qubit gates on 
distant sites, even if it permits complete controllability. 

The specific problem addressed here is the following 
one: we study controllability by the alternate application 
of two Hamiltonians. One of the Hamiltonian describes 
a nearest-neighbour interaction defined by some graph. 
The other Hamiltonian is a projector given by the char- 
acteristic vector of a subset of vertices. This describes 
an interaction between every two spins associated to the 
elements of the subset. 

Our main result is to characterize a large family of 
graphs that give a pair of Hamiltonians implementing any 
quantum dynamics, thereby rendering a system control- 
lable (Section II). The result can be seen as an analogue 
of the Burgarth-Giovannetti infectivity criterion [Tl| in 
this setting involving two Hamiltonians. A comparison 
with the infectivity criterion, possibly with the use of the 
notion of zero-forcing [l2j , remains an open question. It 
is also an open problem to determine the necessity of our 
condition. 

The proofs follow easily from facts of Lie theory and 
algebraic properties of graphs. The members of our fam- 
ily present a particularly disordered cycle structure [l3| . 
Specifically we require that the number of cycles of a 
certain length, starting from different vertices, cannot be 
written as a sum of the numbers of smaller cycles. We 
will show that this is a property responsible for control- 
lability, when the Hamiltonian is the adjacency matrix 
of the graph. Indeed, we will use powers of the adja- 
cency matrix. These encode the cycles of a graph (see 
the definition of a walk matrix below). In analogy with 
a known result in quantum control theory of spin sys- 
tems, the path graph turns out to be the arguably most 
simple example The setting is directly equivalent 

to a single excitation evolving on an XY spin chain with 
constant couplings (here XY means XX + YY, a la Bose 
@). The path is a connected graph with minimal number 
of edges, therefore it corresponds to a very sparse Hamil- 
tonian. This is a fact to take into account, since sparse 
Hamiltonians can be simulated efficiently in a quantum 



2 



computer 14 1. 

In general, focusing only on the dynamics restricted 
to the n-dimensional subspace, the physical device for 
implementation consists of any machine realizing QWs 
(e.g., an optical waveguide lattice [Hj]). While we do 
not modify the circuitry for different tasks, an external 
clock is necessary, since we need to know the time of 
application of each Hamiltonian, even if the resulting op- 
eration is just a phase factor. We will give evidence that 
our property is almost sure (Section II). This is paral- 
lel to the fact that almost every (generic) Hamiltonian 
gives sufficient control. For many types of graphs, the 
construction of an infinite family with a typical property 
may not be straightforward (e.g., expanders, small diam- 
eter graphs, etc.). We will present a method to construct 
infinite families of graphs (without special constraints) 
that satisfy our required property (Section III). 

The results of the paper may translate into valuable 
information in the perspective of designing schemes for 
scalable quantum computation via local control (e.g., 
help in the selection of the systems, the engineering of 
control functions, etc.). (See [3|, [5| for extended treat- 
ments of this topic.) From a wider angle, the results con- 
sist of a step towards a better and more general under- 
standing quantum evolution on networks. Additionally, 
we introduce concepts that propose an interface between 
control theory and graph theory. 

The paper has four further sections: Section II contains 
the general result; Examples are in Section III; we draw 
some brief conclusions and state open problems in Section 
IV. 



II. CONTROLLABILITY 

Let X — (V(X), E(X)) be a (simple) graph with a set 
of n vertices V(X) and a set of edges E(X) C V(X) x 
V(X) — {{i,i} '■ i G V(X)}. The adjacency matrix of 
X, denoted by A(X), is an n x n matrix with A(X)ij 
1 if {i,j} e E(X) and A(X) hJ = 0, otherwise. The 
walk matrix of a graph contains information about its 
cycle structure [16]. Let X be a graph on n vertices 
and let z G R". We define and denote by W Z (X) = 
(z A(X)z ... A(X) n ~ 1 z ) an n x n matrix with entries 
in Z-° associated to X. When z is the characteristic 
vector of some set S C V(X), the matrix W Z (X) is called 
a walk matrix of X with respect to S. In this case, we 
may write Ws(X) instead of W Z (X). Let X be a graph 
on n vertices and let z £ M™. The pair (X, z) is said 
to be controllable if the matrix W Z (X) is invertible (i.e., 
det(W z (X)) ^ 0). When z is the characteristic vector 
of some set S C V(X), we may write (X,S) instead of 
(X, z). The definition of a controllable graph arises as a 
special case of a controllable pair. Let X be a graph and 
let 1 be the all-ones vector. This is also the characteristic 
vector of V(X). The graph X is said be controllable if 
(X, 1) (or, equivalently, (X, V(X))) is controllable. 

Let us recall that a walk of length I in a graph X is a 



sequence of vertices 1, 2, 1,1 + 1, such that {i, i + 1} £ 
E(X), for every 1 < i < I. The ij-th entry of the 
walk matrix, [W x (X)\ hj = Y%=i ( x )i counts the 
number of all walks of length I — 1 from vertex i. Let 
'■— \{j '■ {hi} £ E(X)}\ be the degree of a vertex i. 
A graph X is regular if d(i) is constant over V(X). No- 
tice that a controllable graph can not be regular. In fact, 
the walk matrix of a regular graphs has rank 1, because 
1 is one of its eigenvectors. One can verify by exhaustive 
search that there are no controllable graphs on n < 5 
vertices. Fig. 1 below contains drawings of all connected 
(non-isomorphic) controllable graphs on six vertices. Nu- 
merics show that the ratio [number of graphs] / [number of 
controllable graphs] decreases with n (see caption of Fig. 
1 for small examples). We expect that asymptotically al- 
most surely every graph is controllable, also considering 
that the authomorphisms fixing the vertices of a control- 
lable graphs are trivial. 




FIG. 1: Drawings of all connected non-isomorphic control- 
lable graphs on six vertices. The vector (d(i) : i G V(X)) 
is the degree sequence of X. The degree sequences of these 
graphs are particularly irregular. First line, from left to right: 
(1,2,2,2,3,4), (1,1,2,3,3,4), (1,1,2,2,3,3), and (1,2,2,3,4,4). Sec- 
ond line, from left to right: (1,2,2,3,3,3), (2,2,2,3,3,4), 
(1,2,3,3,3,4), and (2,2,3,3,4,4). There are exactly 8,85,2275 
(connected) controllable graphs on 6, 7, and 8 vertices, respec- 
tively. 



A (continuous-time) quantum walk on a graph X, 
starting from a state |?/>o) £ C™, is the process induced 
by the rule U M (x){t)\i>o) ^ \i>t), where U M{x) (t) := 
e -iM{x)t e |+) anc j M(X) is a symmetric matrix with 
nonzero entries corresponding to the edges of X (e.g., ad- 
jacency matrix, combinatorial Laplacian, etc.). A proba- 
bility distribution supported by V(X) is obtained by per- 
forming a projective measurement on the state \tpt)- The 
matrix M(X) can also be seen as governing the dynam- 
ics of a system of spin-half particles restricted to a single 
excitation subspace. The dimension of such a subspace 
is in fact n. Here we work with adjacency matrices only, 
but the results described are valid for any symmetric ma- 
trix. Studies of perfect state transfer and entanglement 
transfer in spin systems are often carried on with respect 
to this restriction [TtJ ■ QWs and their discrete analogues 
(e.g., coined quantum walks, scalar quantum walks, etc.) 
have found a number of algorithmic applications. The re- 
views p} give a detailed perspective on this and related 
topics. 
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When choosing Hamiltonians of the form of M(X), the 
question to ask about controllability is the following one: 
can we obtain any quantum dynamics on an n-level sys- 
tem by performing repeated applications of QWs? More- 
over, how much can we limit our resources (e.g., number 
of non-null interactions, number of different Hamiltoni- 
ans, etc.)? In particular, can we use just a single QW 
(i.e., a fixed graph) plus an extra operation acting on a 
subspace of a relatively small dimension? The latter one 
is linked to core questions in quantum control theory, 
where we are interested in driving global dynamics by 
directly modifying only a limited portion of the system 
under a parsimony criterion. In the quantum mechanical 
set-up, controllability occurs together with the ability of 
constructing with reasonable accuracy any unitary ma- 
trix of the appropriate dimension (see Chapter 3 of [2]). 
The corresponding property is expressed if we guarantee 
density in U (n) of the group of unitaries realized as se- 
quences of QWs. The next technical lemma describes a 
relation between controllable pairs and Lie algebras. 

Lemma 1 Let X be a graph and let z be the character- 
istic vector of a set S C V(X). Let us define the sym- 
metric (0, 1) -matrix L — zz T . If (X, S) is a controllable 
pair then the real Lie algebra generated by the matrices 
A(X) and L is Mat„ xn (R), the algebra of all n x n real 
matrices. The real Lie algebra generated by iA(X) and 
iL is the vector space of all skew-Hermitian matrices. 

Proof. We prove by induction on k that the Lie algebra 
generated by A = A(X) and L contains the matrices 
A k ~ % LA l , with i = 0, . . . , k. The first claim in the lemma 
will follow at once from this. We note that there are 
integers c r such that LA r L = c r L. If our Lie algebra 
contains the matrices A k ~ % LA % , then it contains the Lie 
products A k+1 ~ l LA l - A k ~ l LA l+1 , with i = 0,...,k, 
and the partial sums A k+1 ~ l LA l - LA k+1 , for all i. In 
particular, it contains A k+1 L — LA k ^ 1 ~ % and therefore 
also 

LA k+1 L - L 2 A k+1 - A k+1 L 2 + LA k+1 L 
= 2c k+1 L - C() (A k+1 L + LA k+1 ). 

for all i. From this, it follows that it contains LA k+1 , 
and therefore all the monomials A k+1 ~ l LA 1 . Let us now 
consider the second claim of the lemma. We say a matrix 
is a commutator of degree r + 1, if it can be written as 
AX — XA or LX — XL for some commutator of degree r, 
where the commutators of degree zero are the matrices in 
the span of A and L. Since A and L are symmetric, we see 
that all commutators of even weight are symmetric and 
all commutators of odd weight are skew-symmetric. The 
intersection of the space of symmetric matrices with the 
space of skew-symmetric matrices is the zero subspace, 
from which we deduce that the even- weight commutators 
span the space of real symmetric matrices and the odd- 
weight commutators span the space of skew-symmetric 
matrices. This implies that the even-weight commuta- 
tors in iA and iL span the space of skew symmetric ma- 
trices, with dimension (n 2 — n)/2, while the odd-weight 



commutators span a complementary space of dimension 
(n 2 + n)/2. This proves our second claim. ■ 

We remark that it is not hard to show that the matrices 

(o o) and (! o) 

generate Ma^2x2(K), but the real Lie algebra they gen- 
erate is s£(2,K) rather than q£(2,M.). Lemma [T] is rem- 
iniscent of the Lie Algebra Rank Condition in quantum 
control theory. The following result gives a sufficient con- 
dition to render a system controllable by a QW: 

Theorem 2 Let X be a graph and let z be the character- 
istic vector of a set S C V(X). If (X, S) is a controllable 
pair then the unitary matrices C^4(X)( S ) = e~ lA ( x ) s and 
Uh(t) = e~ lLt , s,t £ R-°, generate a dense subgroup of 
the unitary group U(n), n>2. 

Proof. Let G be the closed subgroup generated by the 
given elements. Then it is a Lie subgroup of U (n), and its 
tangent space is the Lie algebra generated by iA(X) and 
iL. Since X is controllable, by Lemma[TJ this Lie algebra 
is the space of all skew-Hermitian matrices, which is the 
tangent space to the unitary group U(n). It follows that 
G = U(n). m 

If S = V(X) then z — 1 and the Hamiltonian L = J, 
the all-ones matrix. For its unitary, we have Uj(t)k,i = 
I ( n + e -int _ 1) if k = i and Uj(t) k ,i = i (e- nit - l), 
otherwise. In fact J/m^) is a polynomial in M with 
degree at most the degree of the minimal polynomial 
of M. Among adjacency matrices, this is the fullest 
possible Hamiltonian. Its implementation has been dis- 
cussed in several works QJ]. The matrix Uj(t) is essen- 
tially the same as the Grover operator used in quantum 
search algorithms (l9| . Turning our attention to differ- 
ent characteristic vectors, we can prove a similar result 
concerning the path graph. We denote by P n the path 
of length n — 1, i.e., the graph on n vertices {1,2, ...,n} 
and edges {1, 2}, {2, 3}, {n — 1, n} (1 and n are called 
end-vertices). Weighted paths are often used to model 
ID spin chains. A control criterion concerning the global 
space of ID spin chains have been isolated in |5|. Con- 
trollability of these systems and scalable quantum com- 
putation has also been discussed by in Q. 

Corollary 3 Let P n be the path on n vertices. The uni- 
tary matrix UA{p n )(s) — e~~ lA ( p ™) s together with the diag- 
onal unitary matrix U eie r(t) — e~ leiCl 1 — [e~ lt , 1, 1], 

s,t £ M-°, where e\ — (1, 0, 0) T , generate a dense 
subgroup of the unitary group U(n), with n > 2. 

Proof. We need to prove that W(P n ) = 
( ei A(P n )a . . . A(P n ) n ~ 1 ei ) is invertible. Observe 
that the first entry of the vector A l (P n )ei is a Cata- 
lan number C;/ 2 = jjpi (1/2) ^ ' ^ s even an d zero, other- 
wise; the second entry behaves similarly, with C(; + i)/2 = 
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7^3 ((j+i)/2)' ^ or ' oc ^ - ^ or exam pl e , the first two 
rows of W{P 7 ) are (1, 0, 1, 0, 2, 0, 5) and (0, 1, 0, 2, 0, 5, 0). 
Moreover, the matrix W(P n ) is upper triangular. Since 
Ck >■ "Ylii—oCi, for every k > 3, it follows that the rows 
of W (P n ) are linearly independent. ■ 

III. EXAMPLES 

Given S C V(X), the cone of X relative to 5 is the 
graph X s such that V(X S ) = V{X) U {0}, for a new 
vertex 0, and E{X S ) = E{X) U {{0,i} : i e S}. We 
denote by X\i the graph obtained from X by deleting 
the vertex i and all its incident edges. 

Theorem 4 Given a graph X and a vertex 1 € V^(A) 7 
</ie pair (-Xi, {0}) is controllable if (X, {1}) is. 

Proof. We will show that if u G V(Z) for some graph 
Z, then (Z, {u}) is controllable if and only if the charac- 
teristic polynomials (of the adjacency matrices) <fi(Z,t) 
and 4>{Z\u,t) are coprime. From the properties of X\, 
one can prove that (f>{X\,t) = t(f)(X, t) — <fi(X\l, t). From 
this, we deduce that if 4>(X, t) and <f)(X\l, t) are coprime 
then so are <f>{Xi , t) and 4>{X, t) . Now we derive our char- 
acterization of controllability. Assume n — \V(X)\. Let 
ei be the first vector of the standard basis, and let Eg 
denote the idempotent in the spectral decomposition of 
A = A(X) that corresponds to 9 (see pp. 186-187) 
it follows that 

ffiM = [(tl A)-\i = J> - Oy^Eeer. 

We observe that the number of poles in the rational func- 
tion here is equal to the number of eigenvalues 9 such 
that ef Egei ^ 0; in other words, it is equal to the num- 
ber of eigenvalues 9 such that the projection Ege\ ^ 0. 
Note also that this number is n if and only if <fi(X, t) 
and <f>(X\l,t) are coprime. To complete the argument, 
consider the walk matrix W ei (X). By spectral decom- 
position (again) A r e± = Egei, from which it follows 
that the column space of W Cl (X) lies in the span of the 
nonzero vectors Ege\. Since each projection Eg is a poly- 
nomial in A, we conclude that rk(VF) is equal to the num- 
ber of eigenvalues 9 of X such that Ege\ ^ 0. This proves 
our characterization. ■ 

A method based on the theorem can be used to con- 
struct infinite families of controllable graphs: 

Corollary 5 Let (X, S) be a controllable pair. If Y is 
the graph obtained by joining one end-vertex of the path 
to each vertex in S , then Y is controllable. 

IV. CONCLUSIONS 

We have considered controllability and QWs. As a 
technical tool, we have introduced the combinatorial no- 
tion of a controllable pair. A graph and a subset of its 



vertices form a controllable pair, when the structure of 
the graph exhibit a certain type of disorder. The dis- 
order is expressed by the cycle structure of the graph, 
encoded in the entries of powers of the adjacency ma- 
trix. We have proved that a QW involving a such a pair 
renders a system controllable. By this result, we can in 
principle perform universal quantum computation as an 
alternating sequence of QWs on two graphs, or on the 
same graph, but interspersed with phase factors. Includ- 
ing fault tolerance in this picture would encounter hard 
obstacles because of the sensitivity to phenomena linked 
to decoherence and Anderson localization [22]. An issue 
related to the more abstract aspects is the lack of trans- 
parency when trying to design algorithms with a logic 
that requires operations on specific subsystems. We con- 
clude by stating four problems: Problem 1: Let G be a 
subgroup of U(n) which fixes \ip) € C™ (n > 3) and let 
V € U {n)—G. Then, the group (G, V) {i.e., the subgroup 
of U (n) generated by G and V) is dense in U(n) (see, e.g., 
Lemma 20 in [l|). Is there an analogue to this statement 
for dynamical Lie algebras generated by adjacency ma- 
trices? In particular, let z be the characteristic vector 
of S C V(X) and let P be a permutation matrix corre- 
sponding to an automorphism of X. When Pz = z, it 
follows that PW S (X) = W S (X) and then P = 1 because 
det(Ws(X)) 7^ 0. This means that the authomorphisms 
of X fixing S are trivial if (X, S) is controllable. Is this 
the most general condition for controllability? Problem 2: 
Can we lift the combinatorial criterion for controllability 
introduced in this paper to general criteria for controlla- 
bility of spin systems? Problem 3: Study controllability 
by adjacency matrices, when the time of application of 
each Hamiltonian is constrained. Determining relations 
between quantum control by nearest-neighbour interac- 
tion on graphs and classical simulatability of the asso- 
ciated dynamics is an open problem. Problem 4- What 
can be said about controllability by acting only on a con- 
nected induced subgraph? If the number of vertices is 
constrained, the optimum may be difficult to compute. 
This would be parallel to the Burgarth-Giovannetti cri- 
terion whose optimum is difficult even to approximate 
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